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Abstract 

We consider the low-energy collective excitations at finite temperature of Bose- 
Einstein condensed gases (and liquids as well). A most general model- independent 
effective Lagrangian is written down according to a prescription obtained from the 
breakdown of the global symmetry U{1). To show how the theory predicts easily, we 
derive the momentum and temperature dependence of the damping of excitations by 
means of power counting as an example. 
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1 Introduction 



Bose-Einstein condensation (BEC) was achieved in 1995 in a remarkable series of experi- 
ments on atomic vapors [|I|, and has tremendously boosted both theoretical and experi- 
mental studies (for an up-to-date review, see, e. g., Dalfovo et al. [^). Among others, the 
study of collective excitations in BEC gives us a great test on the finite-temperature many- 
body theory of interacting Bose gases that has been developed over last several decades. It 
is well known that excitations associated with BEC have been well studied in the context of 
liquid ^He [Q. The semi-phenomenological Landau hydrodynamic theory is very successful, 
but it essentially leans on an ad hoc postulated Hamiltonian and the roles of Bose conden- 
sation and broken symmetry are not clear there This theory was justified in a sense 
of microscopic arguments by Feynman 0. On the other hand, beginning with the work of 
London |]^, a microscopic theory of superfluid "^He using field theoretical methods has also 
been extensively investigated through a weakly interacting gas model [Q. The most impor- 
tant results obtained from this model are based on the presence of a Bose condensate, which 
breaks down the global symmetry f/(l) of the system. In fact, the concept of broken sym- 
metry originated in particle physics has been brought into the low temperature phenomena 
of superfluidity and superconductivity first by Anderson [Q, and now it is beheved that it 
controls much of dynamics of many-body systems. 

In the microscopic study of a Bose condensed gas, one usually needs a detailed model 
like that of Bogoliubov, Beliaev , or Gross-Pitaevskii to explain the mechanism for a Bose 
condensate, and as a basis for approximate quantitative calculations, but not to derive the 
most important exact consequences of the spontaneous symmetry breakdown. (Of course, 
such a microscopic model itself is an approximation.) Yet, if the general features are in 
fact model independent consequences of the spontaneous breakdown of the U{1) symmetry, 
why can't we derive them directly from this breakdown by finding the most general effective 
Lagrangian of the system? An answer to this question is provided by a standard technique of 
modern quantum field theory |]lOl , called the effective field theory approach to the symmetry 
breaking. It was applied to superconductivity |TT], where the fundamental properties 
of conventional superconductors such as Meissner effect and flux quantization are derived 
directly from the spontaneous breakdown of electromagnetic gauge invariance without using 
a detailed model like that of Bardeen-Cooper-Schrieffer (BCS) theory. 

This paper aims to show that many important properties of low-energy excitations in a 
generic homogeneous Bose system can be essentially determined as the exact consequence of 
spontaneously broken symmetry without introducing unnecessary approximations. Assum- 
ing the U{1) invariance is spontaneously broken, we shall derive the most general effective 
Lagrangian of the system by performing procedures very much analogous to that of pion- 
nucleon effective field theory []13[ in quantum chromo dynamics. The present paper gives 
an alternative derivation of the Popov effective Lagrangian that we used for calculating the 
damping rate of the excitation [0, |T3[. As we shall see, the whole procedure is incredibly 
simple and straightforward, but also gives a feasible way of predicting some experimental 
observables. An example on the damping rate shall be given to show how the theory works 
easily. The purpose of the present work is twofold. First, it offers a deep way of justifying 
the Landau theory of quantum hydrodynamics and explaining why a weakly interacting gas 
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model can yield some meaningful results for strongly interacting superfluid ^He. Second, 
it shows a simple, powerful effective field theory approach to investigating the excitation 
dynamics of a homogeneous Bose condensed system. 



2 Effective Lagrangian 



We start by considering a simple nonrelativistic many-body problem of spinless bosonic 
particles at finite temperature. In units of h = 1 and ks = 1, its Euclidean action functional 
is given by [|1^, |1^ 

pf} POO 

J[^,^t] = I dr I d^xC 



with the Lagrangian density 
1 



1 1 ■ r .n2 

{x)drip{x) — ^{x)dr4''^ (x) —- — V'?/'^(x) ■ V'?/'(a;) — A i/j^xy^p^x^ 



(2) 



where (3 = 1/T denotes the inverse of temperature. Here, we write x = (x, r) and r 
is the imaginary "time". All fields fulfill the boundary conditions as of ip{'x,T + f3) = 
iIj{'x,t). Obviously, the action (0) is Galilean-invariant, and is invariant under the global 
phase transformations of group U{1), 

ipix) exp{iA)ip{x), 

with A an arbitrary constant. This symmetry is known to be completely broken below some 
critical temperature Tc, i.e., a Bose-Einstein condensate appears with {ip) ^ |jl8|. The 
particle density and current associated with the symmetry are J = (i/2m)(-?/'''"V-?/' — V^V'?/'^), 
J4 = i'^'^'^ = ip, with the current conservation idrP + V ■ J = 0. 

According to our general understanding of spontaneously broken symmetries |T^, any 
system described by an action with symmetry group G, when in a phase in which G is 
spontaneously broken to a subgroup H, will possess a set of Goldstone modes, described 
by fields that transform under G like the coordinates of the coset space G/H. In our case, 
there will be a single Goldstone mode described by a (massless) real scalar field that 
transforms under G = U{1) like the phase A itself. The group U{1) has the multiplication 
rule g{Ai)g{A2) = g{Ai + A2), so under a phase transformation with parameter A the field 
(f){x) will undergo the transformation 



0(x) + A. 



(3) 



In low temperature physics, the Goldstone mode is accompanied with another excitation, 
known as the density (order parameter) fluctuation, which we will see to have nearly zero 
frequency in the long- wavelength limit. Both together form a non-trivial irreducible linear 
representation of the group U{1). To see the theory must involve the Goldstone fields, we 
may write all ordinary complex fields as 
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tjj{x) = \J p{x) exp{i(j){x)) 



(4) 
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where the p is the density field. Under phase transformation, the p is invariant while 
transforms according to the rule @. Now, rewritten in terms of fields and p, the Lagrangian 
(0) becomes 

If the symmetry is broken, we have po = \ {ip)o\'^ 7^ where (■ ■ ■)o indicates the expecta- 
tion value over the ground state. It follows that we may write p(x) = po + cr{x), where the a 
field describes density fiuctuations, known also as collective modes. The effective Lagrangian 
for and a then reads 

£ g = (^^dr(p - — - — — A Po + 0^ , 6 

where we have ignored all total derivatives. For later use, we record the current conservation 
in terms of and a explicitly: 

m^a - ^■|'''° + '^>^^' = 0. (7) 

m 

Here is an important point: to derive Eq. (^) it was not really necessary to start with 
the model Lagrangian (^). Indeed, according to our understanding of spontaneously broken 
symmetry, we did not need to start with any specific theory. A familiar example in particle 
physics is the effective theory of pion-nucleon interaction with SU (2) x SU (2) spontaneously 
broken to SU{2). The important thing is that Eq. @ is invariant under the U{1) trans- 
formation. The general theory of broken symmetries (see, for example, Weinberg [1^) tells 
us that, for symmetry group U{1), the most general form of the effective Lagrangian for 
the density fluctuation fleld and the Goldstone fleld must be constructed solely from the 
ingredients a, Vcr, drcr, V0 and dr(f) together with higher derivatives of these objects. There 
are also two additional rules that the Lagrangian must obey. The first rule is the Galilean 
invariance for a non-relativistic system such that the combination of idr(f)— (V0)^/2m must 
always appear together in the effective Lagrangian and likewise idrCr — (l/m)V0 ■ Va [^]. 



This has been understood through investigations by Takahashi and Greiter et al. |£2 
The second rule is the time-reversal symmetry such that the action is invariant under the 
transformation 

— »• —0, r — > — r. 

This is observed from the action (|l|). This symmetry requires that only the even powers of 
the Galilean invariant [idr(J — (l/m)V0 ■ Vcr] be included in the Lagrangian. For instance, 
the possible lowest power of it is [idT-cr — (l/m)V0 ■ Vcr]^. This is equivalent to have terms 
of [(po + o")V^0/m]^ by making use of the current conservation (|^) [pS] . 

Hence, according to the above prescription the most general f/(l)-invariant action func- 
tional takes the following form 



Ieff[(T, (l)]= dr d xC^Q 

JO J ~oo 
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with 



eff 



cr 



idr<^ — 



(V0)^ 
2m 



^(V0)^ - ^ 
2m^ 2 



2m 



3! 



CnlTL 



-a 



idr 



2m 
2m 



c 



2m " 



2m 3!m 
(V0)2 



a 



2m 



+ 



(9) 



The terms indicated by ■ ■ ■ will contain higher powers and/or derivatives of the a and/or 
fields. Any term of a total derivative has been ignored. The coefficients F^, F^, F^, F„, C2, C3, c' 
and Cg have the dimensions of K^, K, K~^, K~^, K~^, K~^, and K~^, respectively, where 
K represents a typical momentum scale that shall be discussed in the next section. (Here, we 
adopt a normalization such that is dimensionless and a has the dimension of K^.) Similar 
effective Lagrangians appeared in different contexts, but were all based on some microscopic 
model. For instance, Popov [jl9| derived it by means of the power expansion of the pressure 
(or equivalently grand potential) in terms of inhomogeneity for Bose gases, Aitchison et al. 

found it equivalent to a time- dependent non-linear Schrodinger Lagrangian for BCS su- 
perconductors at T = 0, and Demircan et al. implied that it could be obtained from the 
Feynman wave function of superfluids. If — (l/m)V0 and a are identified with the phonon 
velocity field v and the density variation p' of Ref . , respectively, one finds the action (|^) 
corresponds to a Hamiltonian 



-— p + -mp V + 
6\m 2 



(10) 



which is one form of the Landau-Khalatnikov hydrodynamic Hamiltonian 26 . Hence, the 



Landau quantum hydrodynamics is the exact consequence of the breakdown of the U{1) 
symmetry. 



3 Power Counting 

Many-body properties can be studied through propagators (or Green's functions). In the 
following, we write and a fields into a real 2-component scalar 

-C) 

with Greek indices ■ ■ ■ (= 1,2) labeling its components. The propagators are defined 
by the matrix 

A„,(. - = (r{*„(.)*,(=.')}> = /in. #Wd.(xp exp J ffl^. .] 

/[Tlx #(a;)c^o^(a;)] exp / fr[0, a] 
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where T denotes a time-ordered product on r, r'. Consider the quadratic part of the effective 
action 



I^Q = I dr I d X 



2m 2 2m 2m 



quad 

d^xd^x'^\x)V{x,x')^{x'), (11) 



2 

where 

T^i^ _ ( [-^V^ + F^md'^]5\x - x') tdj\x - x') \ , . 

The free propagators are given by the inverse of the matrix T>: 

^{x,x') =V-^{x,x'). (13) 

The calculation of propagators is simplified by transforming to momentum basis via the 
following Fourier transformation 

A(x, x') = X: / rf^fcA(A;)e^''-(— ')~^-(^-^'), (14) 

where Matsubara frequencies = 2t[v / (3 {v = 0, ±1, ±2, ■ ■ ■) and the 4-momentum notation 
k = (k, u,^) is used. We then have 



By finding its inverse matrix, the free propagators (see Fig. |l]) are 



(15) 



with the energy spectrum 



6(k) = - 

m 



F<^k2(F.k2 + c2) 

\ 1 - F;(F<,k2 + C2) 



(17) 



We shall consider the low momentum-energy region such that k = |k| <^ k^ = yc2/-Fo-, in 
which the spectrum reduces to the phonon type 

e(k) ~ ck 



with the phonon velocity c = {1 / m) ^ F^C2 / {1 — F'^C2). Notice that the energy spectrum is 
linear in /c, vanishing as /c — 0. 

In the calculation of Feynman diagrams at finite temperatures in Euclidean field theory, 
one encounters the summation over discrete Matsubara frequencies. A standard technique 
(see, for example, |^) is available to perform such a Matsubara summation. The trick is 
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to use a contour integral in the complex energy plane. Let h{uj) be a function of complex 
variable u, analytical on the line Retu = 0, which decreases faster than as \uj\ oo. 

We then have 

iEM^c..) = -£^/M/^M, (18) 

where 



- 1 ' 

and the C is a contour in complex tu-plane encircling all poles of function h{uj) in a positive 
sense (but those of function f{uj) in a negative sense). 

Now consider a general process involving arbitrary numbers of the Goldstone field 
and the density fluctuation field a. We suppose that their energies and momenta and the 
thermal energy (~ T) are all at most of some order K, which is small compared with ko 
defined above. Even though Lagrangians like (^ are not renormalizable in the usual sense, 



we saw in, for example, the pion-nucleon theory |10| that such Lagrangians can yield finite 
results as long as they contain all possible terms allowed by symmetries, for then there will 
be a counterterm available to cancel every infinity. If we define the renormalized values 
of the constants F^, F^, F^', F^-, C2, C3, C3, ■ ■ ■ in C^q by specifying the values at energies of 
order K, then the integrals in momentum-space Feynman diagrams will be dominated by 
contributions from virtual momenta which are also of order K (because renormalization 
makes them finite, and there is no other possible effective cut-off in the theory). We can 
then develop perturbation theory as a power series expansion in K. 

Each derivative in each interaction vertex contributes one factor of K to the order of 
magnitude of the diagram; internal propagators Aii(/c), Ai2{k) (= — A2i(A;)), and A22{k) 
contribute factors of K~^, and a unit, respectively; and each integration volume df'k 
(= d^kduj) associated with the loops of the diagram contributes a factor of K'^. So a general 
connected diagram make a contribution of order , where 

V = -21^ - /x + E diVi + 4L. (19) 

i 

Here di is the number of derivatives in an interaction of type i, Vi is the number of interaction 
vertices of type i in the diagram, and Jx are the numbers of internal lines and (p-a cross 
lines, respectively, and L is the number of loops. There is a familiar topological relation for 
connected graphs: 

2I^ + I^+E^ = Y,(P,V,, (20) 

i 

where (pi is the number of Goldstone fields in interactions of type i and is the number 
of external lines. Eliminating the quantity J^, the two topological equations above give 



E(rf,-0,)l^. + F^ + 4L. (21) 



The important point here is that the coefficient di — 0j in the first term is always positive 
or zero. Hence, with the numbers of loops and external 0-lines fixed, the leading terms are 
those graphs of di — 0j = 0. The interactions that satisfy this condition are of no derivatives 
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Figure 1: Feynman diagrams of propagators, vertices, and the one-loop selfenergy 11. Notations 
are k = (k, Wj,) and = (</), a) with Greek indices a, /3, • • • = 1, 2. 

of the a field, noticing that the Goldstone fields always appear with derivatives. Next, we 
shall see the power counting (^Tj) provides us a tool that ensures that only finite number of 
terms of £gg are necessary. 

4 Low-energy Excitation Spectrum 

Let us now apply this method to the calculation of higher order corrections to the propagator, 

A'(fc) = A{k) + A{k)U{k)A\k) 

where the matrix n(A;) denotes the self-energy connected graphs. The lowest order correction 
shall be from the graphs of one loop with L = 1 in Eq. (pT]), to which only cubic interactions 
of the Lagrangian (j^) contribute. The following are all cubic interactions that concern us: 



, -^ldr(j){V(l)) , ^{tdr(j)) , 

zm 2 31 dim 



a^ and - 'ja'^^r^ 



which give the following vertices 



6^{k + k' + k" 



E 



.^i k-,k'->k"->k I 



m 



P{k-*k'-^k"^k} 



m 



(22) 

where we define 6'^{k + k' + k") = 6^(k + k' + k")6u+u'+u",o , and J2p{- -} indicates the sum 
over three cyclic permutations. (The vertices and the one-loop diagram of 11 are depicted in 
Fig. [l|.) In the following, we will show how some fundamental low temperature properties of 
low energy excitations can be easily derived without complicated calculations. 

The spectrum of excitations is given by the poles of the exact propagators (Green's 
functions) that are solutions to the following equation 



det A'"^(A;) = det(A-^(A;) - n(A;)) = 0. 



(23) 
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From this expression, the spectrum can be obtained by analytically continuing it to values 
of that are not Matsubara frequencies by doing iu^ ^ u + iri {t] = 0+), but of course this 
shall only be done after the Matsubara frequency sum. From the free propagator expression 
([l6|), the (real) leading term of the spectrum is obviously given by a; = e(k), and we shall 
approximate Keu = e(k) as far as only the lowest corrections are concerned. Thus the 
spectrum is given in the following form 

w = e(k) -Z7(k). (24) 

Now our task is to find the leading order of the imaginary part of the spectrum, called the 
damping rate 7. After some elementary algebras, keeping only leading order in k, Eqs. (p!5|) 
and (H) yield 



-iRe |ni2(k, -'iw + ri) - n2i(k, -iuj + ri)]\ (1 - F^Cj)"' (25) 
with Ili2{k) = —1121 (A;). One can verify that Eq. (p5D can reduce to the specific forms 



previously presented in Refs. [|1^, Rather than go on to really calculate all elements 



of the self-energy matrix, we shall try to quickly derive the temperature and momentum 



dependences of the damping by using the power counting formula (|2T|) . 

At T = 0, the typical momentum scale K is just the momentum A; of a excitation that is 
carried into 11 graphs through external lines. From the expression (0), we find 

nn(A;) ~ k^, Uu{k) = -U2i{k) ~ and U22{k) ~ k\ 

so that Eq. ( [25| ) immediately yields 

7T=o(k) cx k'. (26) 

We therefore have determined 7(k) up to some coefficient that only depends on something 
other than k (e.g., the interaction strength). Expression ( p6D agrees with the well-known 
result first derived by Beliaev in 1958 [^] for a weakly interacting dilute gas model by cal- 
culating the Green's function till the second order approximation. However, our result has 
been derived without assuming a weak coupling, so it is valid for general Bose condensed 
liquids including superfiuid "^He. 

For temperature such that cfc -C T ^ c/cq, we can follow the same arguments as we had 
for T = 0, but some additional examination on diagrams is demanded in order to get the 
correct k and T dependences of 7. Each self-energy element Ua/B carrying a momentum k 
yields a contribution of order K^, which for T 7^ can be decomposed into K" = k^T"~^ 
{I = 0, 1, ■ ■ ■ , z/). For temperature regime concerned here, the fc-dependence of the 11 arises 
only from those vertices that are associated with external (input) momenta after taking off 
the (5-function dependence from the vertex expression ([2^). Hence, fc' can be determined by 
simply counting the power of external (input) momenta from each vertex. In other words, I 
is equal to the number of external (f) lines attached to each H graph. Thus, 



Hn(A;) ~ k^T^ Hi2(A;) = -H2i(A;) ~ kT\ and H22(A;) ~ T^ 
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in comparison with those of T = 0. From Eq. (p5[), we find the damping rate satisfying 

7(k) oc kT\ (27) 
Again, in the context of a weakly interacting dilute Bose gas model, this remarkable result 



was first implied by Mohling and Morita [28| and explicitly obtained by Hohenberg and 



Martin and by Popov [0. However, our result is universally true for a generic system 
of scalar bosons in low temperature and low energy, free of the assumption of weak coupling 
and low density. That implies that the result (l27| ) also holds for liquid ^He. 

The low-energy excitation spectrum (p^ is called phonon spectrum and has an impor- 
tant property: the excitation frequency is equal to zero for zero momentum to all orders in 
interactions, which means that the spectrum does not exhibit an energy gap. This is known 



as the Hugenholtz-Pines theorem |^ that was proved in perturbation field theoretical anal- 
ysis. It has played a crucial role in the early study of superfluidity. Here we gain another 
understanding of this property. It is known from both microscopic theories p2| and 



superfluid hydrodynamics |^ that, if T < T^, the quasiparticle excitations have exactly the 
same spectra as collective excitations (described here by the a field) in the long wavelength 
limit. Further, the effective Lagrangian (^ dictates that the Goldstone mode and the 
collective mode a possess the same pole structure up to all orders in all interactions allowed 
by symmetries. This statement is implied by Eq. (^). For the Goldstone theorem protects 
all kinds of excitations in low temperature and low momentum limit from having energy 
gap, their spectrum must be of phonon type. In this point of view, the Hugenholtz-Pines 
theorem can be simply understood as the equivalent statement of the nonrelativistic version 
of the Goldstone theorem [^, |3^ for the special case of a dilute Bose gas. 



We conclude here that the effective field theory approach to the study of collective exci- 
tations in BEG makes the role of the Bose-Einstein condensate (or broken U{1) symmetry) 
evident, and gives us a model-independent effective action that can immediately predict re- 
sults for experiments. We expect the approach to be very productive and fundamental if 
one can conduct it for trapped alkali vapors that are currently under extensive investigation. 
Yet, the inhomogeneity of the systems due to a trap potential shall alter our results. 
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